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ABSTRACT 

Aims. Extra-solar planet search programs require high-precision velocity measurements. They need to determine how to differentiate 
between radial- velocity variations due to Doppler motion and the noise induced by stellar activity. 

Methods. We monitored the active K2V star HD 189733 and its transiting planetary companion, which has a 2.2-day orbital period. We 
used the high-resolution spectograph SOPHIE mounted on the 1.93-m telescope at the Observatoire de Haute-Provence to obtain 55 
spectra of HD 1 89733 over nearly two months. We refined the HD 1 89733b orbit parameters and placed limits on both the eccentricity 
and long-term velocity gradient. After subtracting the orbital motion of the planet, we compared the variability in spectroscopic 
activity indices with the evolution in the radial-velocity residuals and the shape of spectral lines. 

Results. The radial velocity, the spectral-line profile, and the activity indices measured in He I (5875.62 A), Ho- (6562.81 A), and 
both of the Call H&K lines (3968.47 A and 3933.66 A, respectively) exhibit a periodicity close to the stellar-rotation period and the 
correlations between them are consistent with a spotted stellar surface in rotation. We used these correlations to correct for the radial- 
velocity jitter due to stellar activity. This results in achieving high precision in measuring the orbital parameters, with a semi-amplitude 
K = 200.56 + 0.88 m-s~' and a derived planet mass of M P = 1.13 ± 0.03 M Jup . 

Key words, techniques: radial velocities - stars: planetetary system - stars: individual: HD 189733 - stars: activity 



1. Introduction well constrained, thanks to a series of de tailed spectroscopic 

and photometric monitoring programs dBouchv et al. 2005: 

The precision of radial-velocity (RV) searches for exoplanets B a kos et al. 2006bl IWinnet al. 2007t) . Such monitoring con- 

H . depends on the noise induced by photosphenc luminosity van- smdns the orientation of me system and the poss i b i lity of 

ations due to active regions (e.g. Santos et al. 2002). Until there being another p lanet in the system (ICroll et al. 2007 I : 

now, most active stars have been rejected from major RV sur- |Mdler-Ricci et al^OOSl) . For instance, HST photometry ex- 

veys, but these systems might be interesting objects such as duded with Mgh probability the presence of Earth-sized moons 

HD 189733. Planets around active stars are now found by re- and Saturn . type debris rings around the planet? and excluded a 

search instruments such as CoRoT, and need radial-velocity massiye second planet in outer 2 . , resonance (Pont et al 2007) 
follow-up (Alonso et al. 2008). Monitoring the activity of ac- T _ ,„„„„_ . , , ■ „ • ^, . , 

tive planetary-host stars may provide a means of disentangling 897 p, 1S a chrom ° s P hencall y a <f * ve 8tar - This dwarf 

radial-velocity variations due to Doppler motion from the noise star u 18 as « vanable of .tyP e BY Dra > l e - K stars 

induced by stellar activity, i.e. dark spots in the photospheric re- wl * quasiperiodic photometric variations on timescales rang- 

• ... ing from hours to years. Hebrard & Lecaveher des Etangs (2006) 

urs 1 omn • n 1 1 . u j- detected stellar microvariability in Hipparcos photometric data. 

HD 189733 is a well-known planet-hosting system discov- , , , TT -L , ■, ■ . ■ 

, , n , . , , onrv r. N • . „• u- a u Winn et al. (2006) found that HD 189733 exhibits velocity noise 

ered by Bouchy et al. (2005) in a metalhcity-biased search „ , \ ...... , ,. _i J TT 

c . ... i.j • . r, T . , or photosphenc utter with an amplitude of 12 m s . Henry 

for transiting hot Jupiters with the ELODIE spectrograph, then „ , , • , ■ • 

. , ., ,L , , , ^, . . , tt . & Winn (2008) used the optical photometric variability due to 

mounted on the 1.93m telescope at Observatoire de Haute- , , • , • , r ■ , 

n , nlTT) , T . ■ c a * , • , , t . spots and plages carried into and out of view as the star rotates 

Provence (OHP). It is one of the two brightest transiting sys- r , • , . . , r , 

/\r n i\ a t-x. a » e *u » to determine the rotation period of the star, 

tems (V = 7.7) and thus the parameters of the system are r 

' ' Activity is a signature of stellar magnetism. Moutou et 

* Based on observations collected with the SOPHIE spectrograph al. (2007) completed spectropolarimetric measurements of 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), HD 189733 to derive theoretical models of planet-star magnetic 
France, by the SOPHIE Consortium (program 07A.PNP.CONS) interaction with quantitative constraints. They observed activity 
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signatures that varyied mostly with the rotation period of the star. 
They found that the chromospheric-activity index measured in 
the Ca II H&K lines and Zeeman signatures are consistent with 
a moderate field being present at the stellar surface that is modu- 
lated by the rotation cycle. Shkolnik et al. (2007) monitored the 
Ca II H&K lines chromospheric index and associated the index 
variability with an active region rotating with the star. Moreover, 
they identified an additional flare potentially induced by interac- 
tions with the hot Jupiter. 

In this article, we monitor activity using all the parame- 
ters available from high-precision spectroscopic measurements: 
radial-velocity, bisector of the lines, and chromospheric in- 
dices measured in Hel (5875.62 A), Ha (6562.81 A), and the 
Call H&K lines (3968.47 A and 3933.66 A, respectively). 

Section 2 describes the observations and the data reduction 
process. The refined parameters of the HD 189733b orbital solu- 
tion are discussed in Sect. 3. In Section 4, we examine the radial- 
velocity residuals and the bisector behavior, followed by a study 
of the chromospheric activity in Sect. 5. In the last section, we 
analyze the correlation between parameters and how this analy- 
sis may be used to disentangle radial-velocity variation from the 
noise induced by activity-related phenomena. 

2. Observations and data reduction 

Fifty-five radial-velocity (RV) measurements were obtained be- 
tween 12 July 2007 and 2 3 August 2007 with the high-resolution 
SOPHIE spectrograph (Bo uchv et al. 20061) mounted on the 
1 .93-m telescope at the Observatoire de Haute-Provence, France. 
We gathered about two exposures per night to sample correctly 
both the short orbital period (^ 2.2 days) and the longer rota- 
tion period of the star (^ 12 days). This SOPHIE campaign was 
performed simultaneously with a MOS T photometric campaign, 
discussed shortly hereafter. 

The observations were conducted in the high-resolution 
mode (R-A/6A ~75 000). The simultaneous Thorium-Argon 
(ThAr) lamp was used to achieve higher precision in RV mea- 
surements. The typical exposure time was 6 minutes, sufficiently 
long to reach an average signal-to-noise ratio (SNR) per pixel of 
about 80 at A = 5500 A. RVs were determined using a weighted 
cross-correlation method with a numerical G2 mask: this mask 
was constructed from the Sun spectrum atlas including up to 
3645 lines, following the procedure of Baranne et al. (1996) and 
Pepe et al. (2002) implemented in the automated reduction pack- 
age provided for SOPHIE. We removed the first 8 blue spectral 
orders containing essentially noise and the last 3 orders in the 
red part that distorted the RV calculation due to the presence of 
telluric lines. The effective wavelength range used for the RV 
calculation was between 431 1 A and 6551 A. 

The typical photon-noise error of the measurements was 
about 1.5 m s and the mean uncertainty in the wavelength cal- 
ibration was around 1.0 ms -1 . The main limitation of SOPHIE 
originates in the Cassegrain Fiber Adaper (developed for the pre- 
vious ELODIE spectrograph), which is responsible for the guid- 
ing and centering system in the fiber entrance. Empirically, the 
stable stars followed by the SOPHIE Consortium have an ac- 
curacy around 3 m s so we added a systematic error of 3 m 
s~ 1 . Our average error was then around 3 .5 m s~ 1 , which was the 
quadratical sum of the photon-noise error, the uncertainty in the 
wavelength calibration, and an estimation of the SOPHIE cur- 
rent systematics. 

A perturbation in the apparent radial velocities of the star is 
caused by the planet blocking part of the Doppler-shifted light 
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Fig. 1. Radial velocities of HD 189733 derived from SOPHIE 
spectra collected in July and August 2007 as a function of time. 
The refined Keplerian orbital solution is plotted as a solid line. 
Individual error bars are also plotted. 



from the rotating star. The resulting radial velocity anomaly 
due to the transit of the planet in front of the star is called the 
Rossiter-McLaughlin effect. Five measurements were made dur- 
ing the transit time, calculated with ephemerides from Winn et 
al. (2007). We removed them from the study, plus another mea- 
surement of too low SNR (< 20 at A = 550 nm) due to bad 
weather. 

Background Moon light also polluted part of the exposures, 
shifting the RV measurements by up to 40 m s . Indeed, the sky 
acceptance of the SOPHIE fiber is 3 arc-second and the spec- 
tra are contaminated by sky background. Clouds or full Moon 
on barycentric Julian date BJD-2454000 = 311 days increases 
the level of polluting light. The Moon light shifts the measured 
RVs significantly because the Moon apparent velocity changes 
from +7 km s to -9 km s _1 between 12 July 2007 and 23 
August 2007, whereas the mean RV of the star is around -2.3 
km s _1 . The sky background may be recorded using the second 
fiber feeding the SOPHIE spectrograph but in our case, this fiber 
was used by the ThAr calibration lamp. Consequently, we re- 
moved 16 measurements of RVs and bisectors when there was a 
significant correlation peak with RV ^ Moon apparent velocity 
in the second fiber during an exposure of another star in less than 
two hours of the HD 189733 exposures. These spectra were not 
removed from the activity studies in Sect. 5 because the Moon 
light did not pollute strongly the flux in the lines used as activity 
indicators. Finally, 33 measurements of RVs and bisector veloc- 
ity span were studied. These measurements are available in the 
version of Table 1 at the CDS, which contains in its cols. 1, 2, 
and 3, the time in BJD, the RV, and the error in RV, respectively. 

3. Determination of system parameters 

The Keplerian solution was derived for the 33 RV points selected 
in Sect.2. We fixed the orbital period and the time of perias- 
tron, well constr a ined by previous studies, using tra nsit timing 
dPont et al. 2007t ICroll et al. 20071 IWinn et al. 2007h . We used 
the values estimated by Winn et al. (2007), P = 2.2185733 ± 
0.0000019 days and our periastron time Tq was derived from 
their time of transit center T T = 2453988.8034 + 0.0002 [BJD]. 
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i Table 2. HD 189733 system parameters 
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Fig. 2. Radial velocities of HD 189733 derived from SOPHIE 
spectra collected in July and August 2007 as a function of orbital 
phase. The refined Keplerian orbital solution is plotted as a solid 
line. Individual error bars are also plotted. 



Parameter 


Value 


P mt [days] 


11.953 ± 0.009 


M, [M©] 


0.82 ± 0.03 


r> rj n ..„i 

P p [days] 


Z.Zlqj i 55 [nxeaj 


1 LBJUJ 


z4j JVSo.z4o /D LtixeaJ 




Z.tJ J700.0UJt [11 ACQ J 


v sin i, [km s" 1 ] 


2.97± 0.22 c 


if] 


85.68 ± 0.04 d 


e 


^ 0.008 


V [kms- 1 ] 


-2.2765 ± 0.0017 


K [m s~'] (before correction) 


200.60 ± 2.32 


K [m s~'] (after correction) 


200.56 ± 0.88 


M p [M Jllp ] 


1.13 ±0.03 



a Henry & Winn (2008) 

b Bouchy et al. (2005) 

c Winn et al. (2007) 

d Pont et al. (2007) 



SOPHIE RV measurements were then most accurately fitted 
with a Keplerian function of semi-amplitude K = 200.6 + 2.3 
m s -1 , an eccentricity consistent with zero and a mean radial ve- 
locity V = -2.277 km s ~' (Figs. [Hand [2] Ta b le 2), in agreemen t 
with previous results (Bouchy et al. 2005; IWinn et al. 2 006). 
Our determination of the semi-amplitude is more accurate than 
the value reported by Bouchy et al. (2005) and the value evalu- 
ated from the measurements published by Winn et al. (2006). 

A zero value for the eccentricity is satisfactory because one 
would expect the orbit of the hot Jupiter HD 189733b to be 
nearly circular due to tidal effects. Nevertheless, Knutson et al. 
(2007) discussed the possibility of a weak eccentricity e cos to = 
0.0010 + 0.0002, where a> is the longitude of pericenter for 
HD 189733b, to explain a time delay of 150 + 24 s in the pho- 
tometric observations of the secondary transit. Similarly, Winn 
et al. (2007) determined the limits the orbital eccentricity and 
found that eccentricity might be as high as e « 0.02 if a> « + 

o 

90 . 

We assessed for the detectability of non-zero eccentric- 
ity for the orbit of HD 189733b in the RV data set. We 
used all published RV measurement s, a total of 151 values, 
dBouchv et al. 20051 IWinn et al. 2007b in addition to our own, 
excluding all measurements for data acquired during the tran- 
sit. We fitted the data with a long-term linear-velocity gradient 
v, taking into account that HD 189733 has an M dwarf stellar 
companion (Bakos et al. 2006a) with an orbital period of ~ 3200 
years and an RV amplitude of ~ 2 km s . The trend identified 
between the three datasets was v =1.4 ± 2.3 m s _1 yr . This is 
consistent with the presence of the stellar companion and does 
not require a long-period planetary companion to be explained. 
Leaving e and u> as free parameters, we determined that the ec- 
centricity was lower than 0.008. The long-term velocity gradient 
does not have a significant effect on the determination of the ec- 
centricity. The result agrees with the value of Winn et al. (2007) 
and is still consistent with the Knutson et al. (2007) solution in- 
troduced to explain the secondary transit delay. 

4. RV residuals and bisectors of the lines 

The large "observed-minus-calculated" (O-C) residuals around 
the best-fit Keplerian solution are usually explained by activity- 
induced "jitter" of the star. Our dispersion value, 9.1 m s _1 , is 



smaller than the 12 m s _1 measured by Winn et al. (2007) and the 
15 m s _1 residuals found by Bouchy et al. (2005). These larger 
residuals might be due to the amplitude of the stellar-activity 
jitter decreasing with time (within long-term stellar-activity cy- 
cle) and to the lower accuracy of ELODIE spectrograph used by 
Bouchy et al. (2005). 

The cross-correlation function (CCF), calculated for measur- 
ing the RVs of the star, corresponds to an average of all the spec- 
tral lines. We analyzed the CCF bisector velocity span (V span ), 
a measurement of the shape of the CCF as defined in Queloz 
et al. (2001). Active regions at the stellar surface, such as dark 
spots or bright plages, induce an asymmetry in the emitted light 




Fig. 3. The residuals (O-C) from the Keplerian fit versus bisector 
span of the CCF profile. We observe an anticorrelation, which is 
expected in case of RV variations induced by stellar activity. The 
probability that the 34 observations come from a random sample 
of two uncorrelated variables is lower than 0.1 %. 
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of the star changing periodically as they move in and out of 
view with stellar rotation. They then also induce an asymmetry 
in the line profile that can be measured by the V spa „. An activ- 
ity phenomenon is characterized by an anticorrelation between 
the ob served radial- veloci ty (here, the residuals from the fit) and 
V span dOueloz et al. 200 lh . In Figure. [3] the V spim is anticorre- 
lated with the residuals in agreement with a stellar phenomenon 
as an explanation of the RV variations. V span error bars corre- 
spond to twice the size of the RV error bars. The calculation of 
the V S p an needs two measurements of RV, which provides an ap- 
proximate estimate of its uncertainty. The linear correlation co- 
efficient is -0.58, implying that the null-hypothesis, i.e. that the 
34 observations come from a random sample of two uncorre- 
cted variables, has a probability of lower than 0. 1 %. A linear fit 
to the relation of V span versus residuals (O-C) is given by V span 
= (-0.61 ± 0.15)(O-C)-(12.6 + 1.4). 



5. Spectroscopic indices of stellar activity 

5.1. Call H&K lines 

The Ca II H and K lines are broad, deep lines for solar-type stars. 
Respectively centered on 3968.47 A and 3933.66 A, (Figgl, the 
emission reversal in the core of the Ca II H&K resonant lines, 
resulting from non-radiative heating of the chromosphere, is in- 
dicative of spots and plages, i.e. active regions. When study- 
ing the Ca II H&K lines in SOPHIE spectra, one must account 
for the fact that they are located in a spectral region that suf- 
fers from contamination from the simultaneously recorded ThAr 
spectrum. The star spectrum is recorded in fiber A, while the 
ThAr lamp calibration spectrum is recorded simultaneously in 
fiber B. The calibration lamp produces a background of scat- 
tered light all over the CCD at the level of 0.1% of the stellar 
continuum. In order to use the signature of flux emitted in the 
Ca II lines, we substract this background level from each stel- 
lar spectrum. W e define our Ca II H+K index follo wing the Mt. 
Wilson S index (I Wilson 19681; iBaliunas et al. 19951) : 



Index 



II ■ K 
B+V 



(1) 




a 0.001 

0.0005 

3900 3920 3940 3960 3980 4000 

Wavelength [ A] 

Fig. 4. Ho- (top left) and He I (top right) and emission reversal in 
the Call H and K lines of the star HD 189733 (bottom) in the 
SOPHIE spectrograph typical spectrum. 




X 0.48 
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C 0.475 
- 0.47 
0.465 



300 310 320 330 

BJD - 2454000 

Fig. 5. Evolution of our Ha (top), Ca II (middle) and He I (bot- 
tom) activity indices as a function of time. Error bars are also 
plotted. The red dashed lines indicate roughly the periods of 
minimal and maximal stellar activity based on the Ho- index and 
separated by half the rotation period of the star. The blue filled 
points come from spectra rejected in the RVs and V span studies 
(see text). 



where H and K are the flux measured in a 1 .09 A wide window 
centered on each line of the Ca II doublet, and B and V estimate 
the continuum on both sides with the flux measured in 20 A wide 
windows respectively centered on 3900 A and 4000 A. We de- 
termined the errors of the Ca II index by taking into account the 
CCD readout and photon shot noise. We were unable to estimate 
the error due to the background substraction because it depends 
on too many parameters: SNR, exposure time, elevation of the 
target, and the weather. So, this index is mainly dominated by 
systematics. 

We removed from this study all spectra for which the Ca II 
spectral range had insufficient flux due chiefly to the atmospheric 
refraction that dispersed the blue part of the spectra more signifi- 
cantly. We plot the calculated index for the remaining 38 spectra 
as a function of time in Fig. |5](middle panel) and observe a de- 
crease in the activity level around BJD - 2454000 = 304 days. 



5.2. Ha line 

The Ha line at A = 6562.808 A (Fig. [p is similarly sen- 
sitive to the mean chromospheric activity flKiirster et al. 2 003: 
iBonfils et al. 20071) . Ha is less affected by the flux of the ThAr 
calibration lamp than Call lines. We defined our index as de- 
scribed in Bonfils et al. (2007): 



Index ■ 



'Ha 



F\ +F 2 



(2) 



with F[]a sampling the Ha line and F\ and F 2 the continuum 
on both sides of the line. The F n a interval is 0.678 A wide and 
centered on 6562.808 A, while F\ and F 2 are integrated over 
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0.2 0.4 0.6 0.8 

Planet orbital phase 

Fig. 6. Nightly averaged Ha index as a function of rotation 
phase fitted by a sinusoidal function (top panel). Residuals from 
the individual Ha index after this correction are plotted as a 
function of the planet's orbital phase (bottom panel). In spite 
of a large dispersion due to the short-term activity of the star, no 
specific correlation with the planet motion is observed. 



[6545.495 - 6556.245] A and [6575.934 - 6584.684] A, respec- 
tively. The errors bars are also calculated in a similar way. 

We used all 55 available spectra. Figure [5] shows the varia- 
tions in the Ha index with time; it is characterised by a lower 
noise level than the Ca II index, as expected, and by two dips in 
the activity level around BJD - 2454000 = 303 and 316 days. 

5.3. He I line 

The absorption line of the He I D3 triplet at X = 5875.62 A 
(Fig. [4]) is a signature of a non-radiative heating processes. The 
He I is known to be a tracer of the solar plages (Landman 198 1). 
Since it is sensitive to the heating of the chromosphere (Saar 
et al. 1997), this absorption line is also a diagnostic of stellar 
magnetic activity. We define the He I index in a similar way to 
that for Ha : 



Index — 



'Hel 



F x + F 2 



(3) 



with Fnei sampling the Hel line, and F\ and F2 representing 
the continuum on both sides of the line. Our Fn e i is centered on 
A = 5875.62 A and is 0.4 A wide, while F\ and F2 are integrated 
over [5874.0 - 5875.0] A and [5878.5 - 5879.5] A, respectively. 
The Hel line is located at the borders of two spectral orders 
in our SOPHIE spectra and is far noisier than the Ha line. 
Moreover, in this region of the spectra, telluric lines increase 
the noise level. Nevertheless, two decreases of the index around 
BJD - 2454000 = 302 and 315 days are evident in Fig. [5] and 
perhaps a partial third one around 328 days. The behavior of 
this index reinforces the active region interpretation. 

Shkolnik et al. (2007) identified the presence of an additional 
variability or flaring in the Call cores of HD 189733, potentially 



induced by interactions with the hot Jupiter. We also used the Ha 
activity index to study the chromospheric variability as a func- 
tion of the planet's orbit, since it is less affected by noise than the 
Ca II index in our data. Removal of the main modulation induced 
by stellar rotation was completed by estimating a sinusoidal fit 
of the nightly averaged Ha indices. We did not detect significant 
modulation in the residual activity index with the orbital phase, 
as shown in Fig. [6] This result differs from the tentative conclu- 
sion by Shkolnik et al. (2007), which was based on a 4-night 
data set. It is possible that Ca II and Ha indices are not strongly 
correlated in this star, which may explain the different behavior 
that we observe. Alternatively, it could be another demonstra- 
tion of the on/off nature of star-planet interactions (as argued by 
Shkolnik et al. 2007). 

6. Discussion 

6. 1 . Periodicity of the parameters 

The period derived from an active region is not expected to 
be exactly the rotation period of the star. As for the Sun, 
HD 189733 exhibits differential rotation, with the poles rotating 
more slowly than the equator (Moutou et al. 2007, Fares et 
al., in prep.). We had insuficient data to determine exactly the 
period of our indices of activity. So, we were unable to estimate 
the latitude of the main active region. The rotation period of 
the star has been most reliably determined by Henry & Winn 
(2008). They measured P rot = 1 1.953 + 0.009 days by analysing 
the quasiperiodic photometric variability observed during about 
one year. In Fig. |7l both the residual RV jitter and the bisector 
span appear to exhibit periodic variations on a timescale to 
the rotation period of the star of 12 days. This periodicity was 
consistent with consideration of the parameters as a function 
of stellar rotational phase. We excluded the data points after 
BJD - 2454000 = 320 days because data points in this time 
interval were not in phase with the previous ones. We may 
infer that the main active region moved at this period. This 
is in agreement with the stellar classification for which spots 
are known to evolve quickly. In the same Fig. [7] the lowest 
levels of activity in the Ha line are separated by about 12 days. 
The 12-day variation exhibited by the Ha index becomes more 
obvious when plotted as a function of the rotational phase. We 
excluded the data points after BJD - 2454000 = 320 days that 
had clearly a higher level of activity than for the days before. 
The spectroscopic line profile, the indices of chromospheric 
emission, and the radial-velocities (O-C) therefore exhibit 
periodic variability around the rotation period of the star. 



6.2. Correlations of the parameters 

By searching for the variations in the different parameters to- 
gether in Fig. [7] we can conclude that they are consistent with 
one main active region on the stellar surface, as a first approx- 
imation for multiple active regions. Henry & Winn (2008) esti- 

o 

mated a lower limit of i* ^ 54 for the inclination of the stellar 
equator with respect to the sky plane. This is consistent with 
the fact that the timescale of the maximum-minimum activity is 
around half a rotation period. 

We selected the Ha index as an indicator of activity vari- 
ability since it was the most sensitive index. The activity vari- 
ations exhibit roughly a quarter-phase shift with respect to the 
radial-velocities and line profile variati ons. This fact is in agree- 
ment with an observed stellar activity (Oue loz et al. 20 01) and a 
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Fig. 7. Left panels: Residuals from the Keplerian fit (top), bisector span of the CCF profile (second panel), activity index measured 
in the Ha line (third panel) and MOST photometric data (bottom) as a function of time. Right panels: Same parameters as a function 
of the stellar rotational phase, excluding data at BJD-2454000 > 320 days. In the right panel, a linear trend is substracted from the 
MOST data (see text). The To chosen for the rotational phase is the transit epoch of Winn et al. (2007): Tj = 2453988.8034 days and 
P ro , = 11 .953d. The unit of the MOST data is instrumental magnitude. V spa „ error bars correspond to two times the RV error bars. 
The red dashed lines are based on the minimum activity index in Ha line separated by 6 days, approximately half of the rotation 
period of the star. The plus (minus) signs indicate when the star is at its maximum (minimum) level of activity. Periodic variations 
can be seen for all parameters. (O-C) and V span are antic orrelated (Fig. [3) and the activity index variation is quater-phase shifted 
compared to the RV as expected in case of spot-induced RV variations (Fig. [8j. The photometry is also in line with the activity 
interpretation as the star becomes darker (more spots visible) during the phases of higher activity. 



model of a dark spot rotating with the stellar surface and moving 
in and out of view. In Figure|8] we also plot the residuals from the 
Keplerian fit (O-C) as a function of the chromospheric activity 
index. We observe a loop pattern as expected for an active region, 
which features a dark spot rotating on the stellar surface as al- 



ready described by Bonfils et al. (2007). The quarter-phase shift 
and the loop pattern results from the same correlation. When the 
activity level is at its maximum and the spot is on the central 
meridian of the star, it occupies symmetrically the approaching 
and the receding halves of the stellar disk, and we do not ex- 
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Fig. 8. Radial velocities of HD189733, corrected for the signa- 
ture of its known planetary companion as a function of the Ho- 
spectral index. The loop pattern is a signature of a stellar surface 
with a main active region in rotation. 

pect the RV to be perturbed. Similarly, when the spot is out of 
sight, the activity level is at a minimum and the impact on the RV 
should be negligible, in addition to the impact on the line shape. 
On the other hand, when the spot appears or disappears on the 
visible stellar disk, it lies on the rotationally blue-shifted or red- 
shifted half part of the stellar surface and induces an asymmetry 
in the velocity distribution of the emitted flux. Thus, when the 
dark spot appears (disappears), the activity index is at its middle 
level and the RV effect is minimum (maximum), while the V span 
is anticorrelated with the RV. 

6.3. MOST photometric data 

Our observations of HD 189733 can be compared with those 
taken simultaneously with the MOST microsatellite. MOST 
(Walker et al. 2003; Matthews et al. 2004) is an optical mi- 
crosatellite that is designed to observe nearby, relatively bright 
stars through a single broadband (350-700 nm) filter using a 
15/17.3 cm Rumak-Maksutov telescope. MOST observed HD 
189733 for over ~31 days in 2007. The 2007 observations, 
reduction and resulting analysis is discussed in Rowe et al. 
(in prep.). For comparison with the SOPHIE observations, we 
present the MOST observations of HD 189733 in 2007 in Fig. □ 
(bottom panels). The MOST lightcurve has the transits of the 
known planet removed, the data is binned to the orbital period 
of the satellite (P ~ 101 .43 minutes). A long-term linear trend of 
5.7xl0~ 4 mag day 1 is removed in the right panel of the Fig. [7] 
in order to emphazise the correlations with the other parameters. 
The comparison stars do not exhibit this drift, suggesting that the 
drift is real. Since the timescale for the drift appears to be much 
longer than the length of the observations, one can only specu- 
late to its cause. As can be seen in the bottom panel of Fig. [7] HD 
189733 displayed ~1% variations in flux that were most likely 
due to starspot modulation over -2.5 rotation cycles of the star 
in optical light during the MOST observations. In 2006, during 
the same period, MOST observations of HD 189733 displayed 
~3% variations in flux (Croll et al. 2007). The difference in the 
amplitude of the variation may be due to a decrease in amplitude 
of the stellar activity jitter with time due to a long-term activ- 
ity cycle of the star. This interpretation is only speculation but 
would agree with the decrease in the (O-C) residuals around the 
fit in RV observed by the successive RVs campaigns in 2005, 
2006, and 2007 discussed in Sect. 4. 

MOST photometric data exhibit variations in the stellar flux 
with a periodicity of about 1 1 days. When the MOST data were 
acquired simultaneously with the SOPHIE data, the maximum 
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Fig. 9. Profile of the difference between a mean CCF computed 
at maximum activity level and a mean CCF of all the data avail- 
able (solid line) and the difference between a mean CCF com- 
puted at minimum activity level and a mean CCF of all the data 
available (dashed line). This illustrates the behavior of the spec- 
tral lines with stellar activity. 



flux agrees with the minimum of activity, showing that the vari- 
ation in emitted flux of the star is dominated by dark spots. 
Between BJD-2454000 = 318 and 329 days, the photometric 
curve seems to present a new main active region. The mini- 
mum flux at BJD-2454000 = 319.5 days comes just 10 days after 
the previous minimum and is also deeper than the previous one. 
Moreover, a third shallow minimum occurs at BJD-2454000 = 
324 days, and the maximum flux at BJD-2454000 = 329 days 
comes later (12 days) than the previous maximum. This corre- 
sponds to the fact that the radial-velocities residuals, W span , and 
Ha index data after BJD-2454000 = 328 days are not in phase 
with their previous values before BJD-2454000 = 320 days. 

6.4. Profile of the CCF 

Using the fact that the CCF is equivalent to an average of all 
spectral lines, we compared the profile of the CCF when the star 
was both at its maximum level of activity and its lowest activ- 
ity level. The abscissa of each CCF was first corrected for the 
orbit of the planet. All the available spectra unaffected by the 
background Moon light were used to compute the mean CCF. 
On the other hand, we summed selected spectra to compute a 
mean CCF at maximum activity level (MaxCCF) and a mean 
CCF at minimum activity level (MinCCF). The mean CCF was 
subtracted from MaxCCF and MinCCF. The resulting profiles 
are plotted in Fig. [9] When the star was at its maximum level of 
activity, the difference (MaxCCF-Mean) was positive whereas 
at its minimum level of activity, the difference (MinCCF-Mean) 
was negative. This result could be explained by the presence of 
dark spots when the star was at its maximum level of activity. 
Dark spots have a lower temperature and exhibit different ab- 
sorption lines than those of the star. The spectrum of the dark 
spot may be closer to that of an M-dwarf. When there are dark 
spots on the stellar surface, the lines of the K2V star are shal- 
lower. This result shows that the parameters of the CCF may be 
used to characterize the variation in the stellar activity. 

6.5. Improvement of RV measurements 

In a similar way to Melo et al. (2007), we corrected the radial- 
velocity measurements for the stellar activity and refined the 
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orbital parameters of the planetary system. We used the least- 
square fit between the residuals (O-C) and V span according to 
Fig. |2 as determined in Sect. 4. We added the least-square fit 
between (O-C) and W spa „ to the residuals (O-C) data points. We 
then adjusted the most accurately determined Keplerian fit to the 
corrected radial velocities (new residuals summed to values of 
the previous orbit). As in Sect. 3, we fixed the orbital period, the 
time of periastron, and the eccentricity. We derived lower resid- 
uals from the fit of 3.7 m s _1 , compared with 9.1 m s _1 before 
correction. These lower residuals are in agreement with the esti- 
mated current systematics on SOPHIE measurements discussed 
in Sect. 2. The error bars on the semi-amplitude of the orbit are 
also smaller with K- 200.6 + 0.9 m s _1 , compared with K- 
200.6 + 2.3 m s _1 before correction. Using the mass of the star 
M*= 0.82 + 0.03 M Q given by Bouchy et al. (2005), we found 
a planet mass of Mp- 1.13 + 0.03 Mj up . Our result agrees ex- 
actly with the value of 1.13 + 0.03 Mj up reported by Winn et 
al. (2006), and it is also in agreement with the Bouchy et al. 
(2005) value of 1 . 15 + 0.04 M Jup . The current uncertainty in M P 
is dominated by the error in M*. 

7. Conclusion 

We have observed the transiting planetary system HD 189733 
with the high-resolution SOPHIE spectrograph, with the aim of 
correcting the parameters of the planetary orbit for the effect of 
radial-velocity variations due to activity on the stellar surface. 
We refined the HD 189733b orbital parameters and placed limits 
on both the orbital eccentricity and long-period variation. After 
subtracting the orbital motion of the known planet, the radial- 
velocity residuals were correlated with the spectroscopic activity 
indices and the shape of spectral lines. Chromospheric activity 
tracers and surface photometry were modulated mainly by the 
12-day rotational period of the star. However, the relation that 
we have derived between the bisector span and radial-velocity 
residuals is the most accurate yet and allows a precise correction 
for stellar jitter. The error in the derived semi-amplitude of the 
orbit is reduced by a factor of 2.5 after applying this correction 
law and reduced to close to the instrumental limit, allowing a 
more precise estimate of the companion's mass. 

This result is promising for radial-velocity surveys, which 
cannot continue to avoid targeting active stars. Additionally, 
radial-velocity instruments have to perform follow-up observa- 
tions for photometric planet search programs, which are less sen- 
sitive to the stellar activity for detecting transits. In case the use 
of a single correlation is insufficient, the use of several indices 
allows one to be more confident of accounting for the main ac- 
tive regions and provides stellar-variability diagnostics. In the 
present study, the bisector span was found to be the most reliable 
parameter. In the case of the bright star HD 189733, this method 
of correction is a good compromise. But, if the correction de- 
creases the activity jitter, it also increases the V span photon- 
shot noise in the corrected radial-velocity measurements. This 
method has therefore to be tested further on other systems, to 
estimate the confidence in the corrected RVs and, subsequently, 
to apply similar analyses in planet-search programs. Finally, a 
clearer understanding of the stellar intrinsic variability will have 
a positive impact on identifying and studying activity features 
potentially triggered by star-planet interactions. 
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Table 1. HD 189733 RV measurements 



BJD [day] 


RV [km.- 1 ] 


RV error [km." 1 ] 


2454294.5084 


-2.0865 


0.0033 


2454298.3704 


-2.2320 


0.0037 


2454298.6022 


-2.1180 


0.0034 


2454299.6075 


-2.3931 


0.0034 


2454300.3518 


-2.3702 


0.0034 


2454300.6136 


-2.2123 


0.0035 


2454301.3613 


-2.1437 


0.0034 


2454302.3864 


-2.4464 


0.0035 


2454303.6148 


-2.1555 


0.0035 


2454304.3853 


-2.4705 


0.0035 


2454306.3676 


-2.4206 


0.0038 


2454306.5803 


-2.4638 


0.0036 


2454306.6077 


-2.4644 


0.0033 


2454307.3669 


-2.1593 


0.0034 


2454307.4946 


-2.1057 


0.0033 


2454307.5826 


-2.0795 


0.0033 


2454308.3698 


-2.3264 


0.0033 


2454308.5947 


-2.4465 


0.0033 


2454309.4000 


-2.2859 


0.0034 


2454309.5121 


-2.2181 


0.0033 


2454315.5672 


-2.4759 


0.0038 


2454316.3922 


-2.1139 


0.0034 


2454316.5673 


-2.0794 


0.0034 


2454317.4941 


-2.4306 


0.0034 


2454317.6281 


-2.4652 


0.0034 


2454318.3294 


-2.2216 


0.0036 


2454318.5691 


-2.1097 


0.0034 


2454328.3935 


-2.3662 


0.0033 


2454328.5365 


-2.4255 


0.0042 


2454330.3345 


-2.1936 


0.0037 


2454334.3945 


-2.0675 


0.0034 


2454335.3983 


-2.4625 


0.0035 


2454336.4651 


-2.0804 


0.0034 



